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^L(< Abstract. Polarization observations at A 21 cm and /1 18 cm towards the local Taurus molecular cloud complex were made with 

the Effelsberg 100-m telescope. Highly structured, frequency-dependent polarized emission features were detected. We discuss 
polarization minima with excessive rotation measures located at the boundaries of molecular clouds. These minima get less 
pronounced at the higher frequencies. The multi-frequency polarization data have been successfully modeled by considering 
magneto-ionic Faraday screens at the surface of the molecular clouds. Faraday rotated background emission adds to foreground 
emission towards these screens in a different way than in its surroundings. The physical size of the Faraday screens is of the 
order of 2 pc for 140 pc distance to the Taurus clouds. Intrinsic rotation measures between about - 18 rad nT^ to -30 rad m"^ are 
^ I required to model the observations. Depolarization of the background emission is quite small (compatible with zero), indicating 

a regular magnetic field structure with little turbulence within the Faraday screens. With observational constraints for the thermal 
electron density from Ho- observations of less than 0.8 cm"^ we conclude that the regular magnetic field strength along the line 
of sight exceeds 20 yuG, to account for the observed rotation measure. We discuss some possibilities for the origin of such strong 
■ and well ordered magnetic fields. The modeling also predicts a large-scale, regularly polarized emission in the foreground of 

' the Taurus clouds which is of the order of 0. 18 K at /i 21 cm. This in turn constrains the observed synchrotron emission in total 

, intensity within 140 pc of the Taurus clouds. A lower limit of about 0.24 K, or about 1.7 K/kpc, is required for a perfectly 

ordered magnetic field with intrinsic (~ 75%) percentage polarization. Since this is rather unlikely to be the case, the fraction 
of foreground synchrotron emission is even larger This amount of synchrotron emission is clearly excessive when compared to 
previous estimates of the local synchrotron emissivity. 
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H For large areas of the sky the polarization properties of Galactic . , i • j . . u- u . i . 

d ,. . . „, ■ highly varying polarized structures, which are not correlated to 

■ - - radio emission are not well known in detail. Most information 



cn 
m 



1. Introduction et al. dTWTin^y^ . Gray et al. ( fTWHlIT^^ . Uyamker et al. 

11999), Haverkorn et al. ( 2000), Gaensler et al. fM)l ) revealed 



came from early large-scale surveys of the northern sky caiTied 



structure variations seen in total intensity maps. These findings 

called for a systematic approach and new polarization surveys 

out at Jodrell Bank and Cambridge with 1.4 GHz being the . . . . . 

f ■ were initiated, 
highest frequency used (Bingham [1966). The angular resolu- 
tion of all these surveys was low. Later the Dwingeloo-surveys Here we report on follow-up observations at ^ 18 cm of 
were caiTied out in the frequency range from 408 MHz up to the Effelsberg ^21 cm "Medium Galactic Latitude Survey" 
1415 MHz (Brouw & Spoelstra 1976 1. They have angular reso- (EMLS, Reich et al., in preparation). This survey covers the 
lutions of 0°6 at 141 1 MHz or less at the lower frequencies and northern Galactic plan e in th e range of +20°. The first EMLS 
are absolutely calibrated, but severely under-sampled. Details maps (Uyamker et al. I1999> show an unexpected richness in 
of the polarization structures require fully sampled measure- structure of the polarized emission, which indicates significant 
ments at arc-minute angular resolution. A number of obser- variations in the properties of the magneto-ionic medium. Since 
vations were made during the last decade using synthesis ar- the polarization structures are not associated with total inten- 
rays as well as single-dish telescopes. Polarization observa- sity enhancements or depressions, the interpretation of Faraday 
tions by Junkes et al. (kl987j . Wieringa et al. (1993^, Duncan modulation of an intrinsically polarized background seems the 

most likely explanation for the observations. 

Send offprint requests to: M. Wolleben 

* Based on observations with the Efi'elsberg 100-m telescope oper- The physical parameters of the Faraday rotating regions, 

ated by the Max-Planck-Institut fur Radioastronomie (MPIfR), Bonn, the Faraday screens, are of particular interest as they reflect lo- 

Germany cal conditions of the magneto-ionic interstellar medium, which 
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deviate from the average parameters derived from the analysis 
of rotation measure (RM) from pulsars or extragalactic sources. 
The average density of thermal electrons in the Galactic plane 
is about 0.03 cm"^ and the average of the regular magnetic field 
along the line of sight is about 1 to 2 fiG (Taylor & Cordes 
ITW^ Comez et al. IMJTI Cordes & Lazio ETKBt . However, 
local bubbles with RMs exceeding 100 radm"^ are known to 
exist in the Galaxy (e.g. Gray et al. 119981 Gaensler et al. 
120011 Uyaniker et al. "2003 1. Such large RM values need kpc- 
size structures to originate within the average magneto-ionic 
medium and strongly suggest significant local enhancements 
of the magnetic field strength and/or the thermal electron den- 
sity. The total magnetic field strength was estimated to be about 
5 to 6 juG in the solar vicinity (Strong et al. l2000l l. Even in the 
unlikely case that the regular component along the line of sight 
locally has that value, an excessive size of the Faraday screens 
of some hundred parsec is required. 

The key problem is the unknown distance of the Faraday 
screens and this makes any estimate of their physical parame- 
ters rather uncertain. Distance information can be obtained by 
spectral line observations of either H i, CO or recombination 
lines of neutral hydrogen clouds, molecular clouds or H ii re- 
gions associated with Faraday screens. The Taurus- Auriga re- 
gion seems to be quite suited for such an investigation. The 
large molecular and dust cloud complexes are at a distance of 
about 140 pc, and it is therefore possible to resolve structures 
on pc-scales with arc-min angular resolution. In addition the 
Taurus-Auriga complex is located at medium latitudes well be- 
low the Galactic Plane, so that confusion with the intense emis- 
sion from the Galactic disk is not important. The maps from 
the EMLS revealed a number of polarized enhancements or 
depressions apparently related to molecular material. In order 
to derive physical properties of the associated Faraday screens 
we have complemented the All cm survey data of the Taurus- 
Auriga region by /1 18 cm observations to derive the depolar- 
ization properties and the rotation measures. We modeled the 
polarization data in order to constrain the physical parameters 
of the Faraday screens by taking into account foreground and 
background emission. 

2. Observations and data reduction 

The /1 21 cm observations for the EMLS were recently com- 
pleted and are presently in final reduction stage. The data cov- 
ering the Taurus-Auriga region extending below the Galactic 
plane from I = 150° to / = 190° were extracted and reduced by 
Wolleben J2001t following the survey reduction procedures. In 
brief, total intensities and linear polarization were measured si- 
multaneously with sensitivities of 15 mK for Stokes / and 8 mk 
for Stokes U and Q at an angular resolution of 9'35. The ob- 
serving, reduction and calibration methods of the survey and 
selected maps from different regions along the Galactic plane 
were already published by Uyaniker et al. ( fT998l [T9991 . The 
absolute calibration of the total intensity data was done in the 
usual way by adding the missing large-scale structure from 
the absolutely calibrated Stockert /1 2 1cm survey (Reich & 
Reich"1986). The absolute calibration of the Effelsberg polar- 
ization survey data relied so far on Dwingeloo polarization data 



(Brouw & Spoelstra 19761. However, the Dwingeloo survey is 
rather incomplete in coverage and largely under-sampled, so 
only a few measurements were available for the Taurus-Auriga 
region. We therefore made use of the new polarization survey 
carried out with the DRAO 26-m telescope at 1.41 GHz, which 
consists of fixed declination scans being fully sampled along 
right-ascension. They are separated in declination by 1° to 2°. 
A brief description of the survey methods and preliminary re- 
sults was already given by Wolleben et al. (2004i. The DRAO 
data for Stokes U and Q were compared with the corresponding 
Effelsberg maps both being convolved with a 6°-wide Gaussian 
and the differences were added to the original Effelsberg U and 
Q maps, respectively. 

The A21 cm Taurus- Auriga map in total intensity and po- 
larization was discussed in detail by Wolleben J200H with 
large-scale polarization correction using Dwingeloo data. A 
section of this map with the new large-scale DRAO polariza- 
tion data as described above is displayed in Figs. ^ and b 
showing polarized and total intensities. Total intensities were 
overlayed with the distribution of integrated CO emission as 
observed by Dame et al. j2001t tracing numerous molecular 
gas clouds in this region. 

Follow-up observations of selected fields, where sig- 
nificant polarization structure was seen, were done using 
the Effelsberg 100-m telescope at two frequencies around 
A 18 cm (1660 MHz and 1713 MHz) in September 2001. The 
same receiver as for the All cm EMLS was used, but a dif- 
ferent HF-filter was selected to suppress interference out of 
the A 18 cm band. The observations were done simultaneously 
for both frequencies by utilizing two identical IF-polarimeters. 
Both have an effective bandwidth of 14 MHz. The fields were 
mapped two times in orthogonal scanning directions along 
Galactic latitude and Galactic longitude, in a similar way to the 
All cm EMLS observations. The angular resolution was 7'87 
at 1660 MHz and 7:70 at 1713 MHz, respectively. The maps 
were fully sampled at 3'. 3C286 served as the main calibra- 
tion source with a flux density of 13.3 Jy at 1660 MHz and 
13.1 Jy at 1713 MHz, respectively. The percentage polariza- 
tion was taken as 9.6% at a polarization angle of 33°. In addi- 
tion, 3C138 and 3C48 were used as secondary calibrators. The 
higher frequency maps have been convolved with a Gaussian 
to match the 935 resolution of the All cm EMLS. The rms- 
noise measured were 15 mK (1408 MHz), 16 mK (1660 MHz) 
and 19 mK (1713 MHz) for total intensities and 8 mK, 9 mK, 
and 10 mK for Stokes U and Q, respectively. The same conver- 
sion factor of Tb/S = 2.12 + 0.02 mK/Jy at All cm (Uyaniker 
et al. 19981 applies for all frequencies. 

In order to remove temperature gradients due to varying 
ground and atmospheric radiation from the maps, a linear base- 
line was subtracted from each subscan. This procedure also af- 
fects real sky emission of large spatial extent, this emission 
is usually recovered by absolutely calibrated data. However, 
other than at All cm there exist no absolutely calibrated data at 
/1 18 cm and we calibrated the maps in the following way: The 
temperature spectral index of the Galactic synchrotron emis- 
sion yS' (Tb oc ) was adopted to be /j' x -2.7 in the Taurus 
area following the results of Reich & Reich ( 1988) for the fre- 
quency range between 408 MHz and 1420 MHz. We assumed 
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Fig. 1. a) Section from the All cm "Medium Galactic Latitude Survey" including absolute calibration showing polarized inten- 
sities with polarization vectors overlaid in E-field direction, b) the same for All cm total intensities. These include the 2.7 K 
cosmic microwave background. Overlaid contours show integrated CO emission from the survey by Dame et al. C2001j . The four 
fields observed at /i 18 cm are marked and labeled. 



the same spectral index also for the large-scale (LS) polar- 
ized intensities (JS^^) and calculated an average offset for the 
/1 18 cm maps from that of the All cm map. Rotation measures 
(RMs) across the Taurus region were determined by Spoelstra 
iW71\ using data between 408 MHz and 141 1 MHz. The RMs 
are very low everywhere varying around zero and we assume an 
average /?Mls of rad m"^ to calculate from the extrapolated 
large-scale polarized intensity the offsets for the /1 18 cm U and 
Q maps, respectively. Since the frequency differences are not 
very large we believe that the applied extrapolation closely re- 
sembles absolutely calibrated /1 18 cm polarization data and we 
expect only marginal effects on the small-scale polarized struc- 
ture analysis we are interested in. 

3. Results 

3.1. Polarization maps 

The /i 18 cm maps show - like the All cm map - no remark- 
able structures in total intensities. All maps show smooth dif- 
fuse emission varying mainly with Galactic latitude and a large 
number of unrelated extragalactic sources are superimposed. 



The polarization structures for Fields 1 to 3 (see Fig. ^p) are 
rather similar for all wavelengths with marginal differences in 
structure and polarization angles. Neither of these structures 
show unusual spectral indices of the polarized emission, they 
are not spatially associated with molecular clouds or any other 
known objects and originate in the ISM at an unknown dis- 
tance. For the majority of polarization structures in these fields, 
the present frequency range and separation seems not adequate 
for a reliable analysis of the Faraday effects towards these fea- 
tures. 

In Field 4 there are a number of small-scale polarization 
minima, which spatially coincide with the molecular cloud cov- 
ered in this field. These minima are obviously aligned to the 
boundary of the cloud. They show rather clear differences in 
intensity and thus in spectral index, as well as in the distribu- 
tion of polarization angles between /1 18 cm and All cm. For 
the nine intensity minima we intend to model, we list their 
positions, measured polarization data, spectral index and RM 
in Table I. Spectral index and RM are calculated from the 
1408 MHz and 1660 MHz data. No simple fit based on a /I or A^ 
dependence including all data was possible, reflecting the com- 
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Table 1. Observed values for the center of the identified Faraday screens in Field 4 (see Fig.[TJ)). The observed spectral index 
of the polarized intensity /3^^^ and the observed RMobs were calculated using 1408 MHz and 1660 MHz data. Errors in PI are 
1 1 mK (1408 MHz), 13 mK (1660 MHz), and 14 mK (1713 MHz). 

1408 MHz 1660 MHz 1713 MHz 
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PI 


PA 




PI 
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deg 


deg 


mK 


deg 


mK 


deg 




mK 


deg 








rad m^^ 


a 


170.3 


-9.0 


144 


-29 ±2 


136 


-25 ± 


2 


132 


-23 ± 


2 


-0.4 


±0.9 


-5 ±4 


b 


170.45 


-9.4 


207 


-29 ± 1 


168 


-22 ± 


2 


163 


-21 ± 


2 


-1.3 


±0.7 


-9 ±3 


c 


170.25 


-10.0 


79 


-32 ±3 


104 


-33 ± 


3 


102 


-30 ± 


3 


1.7 ± 


1.3 


1 ±6 


d 


169.75 


-9.7 


97 


-17 + 2 


124 


-28 ± 


2 


118 


-25 ± 


2 


1.5 ± 


1.1 


15 ±5 


e 


169.3 


-9.55 


115 


-5 ±2 


121 


-24 ± 


2 


108 


-23 ± 


3 


0.3 ± 


1 


26 ±5 


f 


171.4 


-8.95 


130 


-18 + 2 


142 


-12± 


2 


147 


-11 ± 


2 


0.5 ± 


0.9 


-8 ±4 


g 


168.15 


-8.6 


75 


-20 ±3 


103 


-30 ± 


3 


90 


-29 ± 


3 


2.0 ± 


1.4 


14 ±6 


h 


174.53 


-10.18 


125 


-77 ±2 


124 


-51 ± 


2 


106 


-55 ± 


3 


0± 1 




-36 ±4 


i 


169.35 


-10.65 


109 


-23 ±2 


117 


-22 ± 


2 


111 


-20 ± 


2 


0.4 ± 


1.1 


-1 ± 5 



plex situation of polarized emission transfer through a Faraday 
screen. Figures |2t and b show the observed polarized intensi- 
ties at 1408 MHz and 1660 MHz with selected E- vectors su- 
perimposed. Figure 121; displays the temperature spectral index 
for the polarized intensity calculated from: 

p, \og{PhlPh) 

log(Vl /V2) 

with the frequencies v and V2 > v\. The spectrum of the large- 
scale emission is close to ffj^^ - -2.7, which was the assump- 
tion made earlier, while the Faraday screens all show flat or in- 
verted spectra with larger intensities at /1 18 cm than at /1 21 cm. 
In Fig. 121 we show the distribution of RM which we have cal- 
culated on the basis of the 1408 MHz and 1660 MHz observa- 
tions. While RMi^s - rad m"^ was assumed for the large-scale 
emission we find clear RM deviations in the Faraday screen di- 
rections. 

4. Discussion 

4.1. Additional data for the analysis 

We searched for additional observations of the Taurus region 
which might be helpful in the interpretation of the polarization 
data. RM variations may result from variations of the thermal 
electron density, which are traced by Ha emission. In Fig.|3t 
we show the Ha emission taken from the "full-sky H-alpha 
map" which Finkbeiner f2003'( compiled from northern sky 
(WHAM and VTSS) and southern sky (SHASSA) data. The 
angular resolution of Finkbeiner's combined map is 6'. No 
structured Ho- emission is visible, which might be correlated 
with the polarized emission features we have measured. The 
quoted Ha sensitivity of 0.52 Rayleigh is taken as a limit for 
associated Ha emission. 

Correlation of polarized emission structures, however, was 
found with integrated CO emission and also with interstel- 
lar absorption. We display a map of the reddening Eb-v in 
Fig. ISJ), which is a reprocessed composite of the DIRBE and 
IRASASSA maps at 100 fim by Schlegel et al. ( IT998t . The 
map of Schlegel et al. is believed to trace the dust density more 
closely than the IRAS 100 fim map. Like the map of integrated 



CO emission these maps show a remarkable correlation with 
the polarization maps. 

4.2. Faraday screens 

The morphology of the observed polarized emission and its 
variation with frequency suggests that strong Faraday rotation 
takes place in regions at the surface of molecular clouds. We 
use a simple model to describe the effect of these Faraday 
screens on the observed distribution of Galactic polarization. 
The screen is placed somewhere inside the synchrotron emit- 
ting region and modulates the polarized emission passing 
through it. According to the total intensity observations, where 
no correlated enhancement is noted, the Faraday screen itself is 
assumed not to radiate and the noise level of the total intensity 
maps provides an upper limit for its emission. 

Galactic synchrotron emission is generated by cosmic -ray 
electrons in the Galactic magnetic field. The position angle 
of synchrotron emission, when corrected for Faraday rotation 
along the line of sight, is in the direction of the magnetic field 
component perpendicular to the line of sight. The two emis- 
sion components in front and behind the Faraday screen are 
supposed to be independent with different intensities and po- 
larization angles, although they might themselves consist of 
multiple layers with different magnetic field orientation and/or 
some internal Faraday dispersion. These emission components 
are assumed to vary on larger scales than the size of the Faraday 
screens and therefore to have the same properties across the 
area modeled. Foreground and background polarization inten- 
sities and polarization angles are described in the following by 
P/foie and PAfoie or P/back and PAback, respectively. 

The Faraday screen may decrease the polarized intensity 
of the background emission by depolarization. In addition it 
changes the background position angle by Faraday rotation. 
Therefore, the observed emission may increase or decrease rel- 
ative to its surroundings, depending on the polarization angle 
orientation of the foreground and the background emission. 

The polarization angle PA of linearly polarized radiation 
is rotated by APA when passing through the magneto-ionic 
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Fig. 2. Field 4: Polarized intensity at 1408 MHz a) and 1660 MHz b) with selected polarization E-vectors superimposed, which 
are shown for every third pixel (separation 12')- The temperature spectral index of the polarized intensity is shown in c), where in 
addition the individual objects being modeled are indicated. The distribution of RM is shown in d). RMs are calculated from the 
polarization data obtained at 1408 MHz and 1660 MHz. Figures|2j; and d show in addition superimposed contours of integrated 
CO emission. 



medium of the Faraday screen and measured at wavelength A. 
The rotation is given by: 

APA =RM-A^, (2) 

with RM, the rotation measure of the Faraday screen: 

RM = Qm B\\nJ . (3) 



Here B|| is the magnetic field strength component along the line 
of sight measured in fiG, is the electron density per cm~^ 
and I is the size of the Faraday screen along the line of sight 
measured in pc. 

In principle a number of effects are able to lower the ob- 
served degree of polarization by depolarization (see Burn I 19661 
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Fig. 3. a) The relevant section from the "full-sky H-alpha map" combined by Finkbeiner ("20031. b) The map of Schlegel et 
al. (1998) showing interstellar reddening Eb-v- Both images show selected contours of the temperature spectral index of the 
polarized emission superimposed. Dashed lines indicate negative spectral indices and solid lines positive ones. 



Sokoloff et al. ll998t . Due to the relatively small spatial extent 
of the Faraday screens discussed here, along with moderate 
electron densities and magnetic field strengths, the following 
mechanisms might potentially contribute: 

Beam depolarization: A variation of the polarization angle 
across the beam reduces the observed degree of polar- 
ization. A rotation measure gradient within the beam can 
cause such an angle variation. ARM - 17 radm^^/beam at 
All cm would reduce the degree of polarization by 10%. 
We do observe such strong gradients towards the Faraday 
screens in our field, and we conclude that beam depolariza- 
tion may not be negligible. Therefore, we include depolar- 
ization in our model. 

Bandwidth depolarization: An interstellar rotation measure 
causes the polarization angle to vary over the frequency 
band of the receiver. A RM of 600 rad m"^ at All cm is 
necessary to decrease the degree of polarization by 10%. 
Such RMs are clearly excessive to those observed here and 
we conclude that bandwidth depolarization is negligible. 

4.3. The model 

In the following a Faraday screen is assumed to be a uniform 
sphere with radius R. Elliptical structures were reduced to cir- 
culars by normalizing the maximum radius at each position an- 
gle of the ellipse to 1. The path length for radiation passing 
through the screen will therefore vary systematically with po- 



sition r. With r as the projection of R on the sky, the path length 
/ is given by l{r) - 1 ^R^ - r^. We also assume a constant elec- 
tron density and a homogeneous magnetic field B inside the 
sphere with a component B\\ along the line of sight. Faraday ro- 
tation is zero for r > R. The modified polarization angle PA^oi 
can be written as 

PA^oiir) = PAback + RMUr) ■ ^ (4) 
with 

RM-Ur) = RMq ■ l(r) ■ (5) 

RMo is the maximum intrinsic rotation measure of the Faraday 
screen at r = 0. 

We assume any intrinsic depolarization DPint to increase 
linearly with the path length l(r). We express int(r) by 

DP^ir) = ^(1 - DPo) + DPo (6) 

with the maximum intrinsic depolarization DPq at r = 0. The 
modified polarized intensity PImod can then be written as: 

Plmod(r) = i^/back ■ DPuJr) . (7) 

In this notation DP - I stands for no depolarization and DP - 
means complete depolarization of the polarized background 
emission. Depolarization is zero outside the screen (r > R). 

Faraday rotation and depolarization will modify the back- 
ground polarization in a systematic way as a function of r (for 
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Fig. 4. The sketch shows the main components of the model. 
Polarized emission from behind the Faraday screen (back- 
ground polarization) is systematically rotated and depolarized 
as a function of path length through the screen. The superposi- 
tion of the modified background and the polarized foreground 
emission leads to a systematic variation of the observed polar- 
ized intensity and angle as seen in the PA-PI plot. 

a schematic representation see Fig. 011. The observed polariza- 
tion is the superposition of the modified background and the 
foreground polarization. In such a case, the polarization prop- 
erties of the combined emission can be computed by separately 
adding Stokes parameters U and Q. Then PI and PA are deter- 
mined from these combined Stokes parameters. Given 



f/mod = PImod ■ sin(2PAmod) 
Gmod = PImod ■ COS(2PAmod), 

the measured polarization becomes 



^/obs('-) 

PAobsir) 



Hon 



+ U^Ur)r + (Gfore + Qmodidr 
iarctan(^^^4^) . 

2 \ Qtoa+Qraod(r) J 



(8) 



(9) 



Polarized intensities and position angles in front of and be- 
hind the Faraday screen depend on the observing frequency. In 
our case we start from absolute polarization measurements at 
All cm and use a spectral index of /S^^ - -2.7 and a /?Mls = 
rad nr^ to adjust the relative /i 18 cm data to absolute values 
as described in Sect. 2. The amount of Faraday rotation by the 
screen is frequency-dependent and affects PA^nod and therefore 
t/mod and 2mod accordingly. Since the beam widths of the ob- 
servations at the three frequencies are rather similar we assume 
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Fig. 5. The plot shows two of the modeled observables for ob- 
ject d at frequencies between 1410 MHz and 1800 MHz. The 
observable rotation measure RMohs (solid), as well as the ob- 
servable spectral index yS^J^ (dashed) are shown as a function of 
the second observing frequency for r = relative to 1408 MHz. 
Vertical lines mark the two /1 18 cm frequencies (1660 MHz 
and 1713 MHz) used in this paper. 



depolarization within the Faraday screen as well as any depo- 
larization DP not to be frequency-dependent. 

Finally the model has four free parameters: P/fore and 
PAfoK (or P/back and PAback, respectively) and the central ro- 
tation measure PMq and depolarization DPq of the Faraday 
screen. 

4.4. Notes on the model 

1 . The PA-PI diagram contains all observed pixels in the field 
of the Faraday screen as outlined in Fig. 5. This also holds 
for the resulting PA-PI relation. 

2. Towards a Faraday screen the observed PMobs differs from 
its intrinsic RMim depending on the amount of foreground 
polarization adding to the rotated background polarization. 
PMobs depends on the two frequencies (vi , V2) used for its 
determination as illustrated in Fig.|5l This also implies no 

dependence of the observed polarization angles in the 
direction of Faraday screens. Figure [S] shows for object d 
how PMobs and spectral index of PI depend on the second 
observing frequency V2, while vi is always 1408 MHz. 



4.5. Application of the model 

For the nine objects identified as minima in the 
/1 21 cm polarized intensities shown in Figs. [5^ and b, 
which also show unusually high 0^^^ as seen from Fig. |2j;, 
we used the following algorithm to apply our Faraday screen 
model and find the best-fit parameters: 

We start by defining the center of the Faraday screen from 
the yS^J^ map for each object and estimate radius and ellipticity 
by eye. This is shown in the left panel of Fig.|6] 

For normalized radii r > 1 no Faraday effects are present 
and the superposition of the foreground and the background 
model component must agree with the observed polarization. 
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Fig. 6. From left to right are shown: The map of fi^^^ in Galactic coordinates with the ellipse marking the size of the Faraday 
screen, the PA-PI plot for 1408 MHz (circles) and 1660 MHz (crosses), the observed spectral index of the polarized intensity 
j3^^^, as well as the observed rotation measure RMob^ versus radius r. Thick lines indicate the modeled values. Each row refers to 
a Faraday screen listed in Tableland labeled from a to i. 



However, in this case Faraday screens may partly overlap and 
influence the off^set polarization. 

With the boundary condition that the model must correctly 
reproduce the observed polarization at r = 1, four parameters 
need to be fitted where we allow variations within a reasonable 
range. These limits are: 50 mK < PI < 500 mK for the fore- 
ground and the background polarized intensity, < DPq < 1 
and -70 rad m ^ < rMq < 70 rad m"^. 



Figure |6l shows the ySj^J^ maps, as well as /3^^^ and RMobs 
versus radius r, and the model fits for the nine Faraday screens 
listed in Tabled All the pixels of a Faraday screen located 
within the ellipse shown in the left-hand panels of Fig.|6lare 
displayed in the PA-PI diagram. We calculated the observed 
spectral index and RM using the 1408 MHz and 1660 MHz 
data. The 1713 MHz data were used to check the model param- 
eters for consistency. 
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Fig. 6. continued 

Table 2. Best-fit parameters of our model fits at 1408 MHz for the Faraday screens listed in Table[l] 



Object 


P/fo,-e [mK] 


/'Afo,, [°] 


Wback [mK] 


/"Aback [°] 


RM(t [radm -] 




a 


186 


3 


179 


-13 


-18 


1 


b 


172 





183 


-10 


-18 


1 


c 


196 


-6 


145 


-11 


-27.5 


1 


d 


200 


-5 


132 


-13 


-28 


0.9 


e 


224 


-2 


120 


-19 


-29.5 


1 


f 


202 


9 


138 


11 


-25.5 


1 


g 


170 


-3 


124 


-13 


-27.5 


0.9 


h 


82 


6 


212 


-15 


-25 


0.9 


i 


172 


-6 


92 


-17 


-22 


1 



In order to estimate how sensitive the model is to system- 
atic effects on the calibration of the large-scale emission of 
the /1 18 cm data, we re-calibrated the data for different spec- 
tral indices and 7?Ms for the large-scale polarized emission 
O^cai = ~2.7 + 0.3 and RM^^x = + 5 radm"^). Subsequent 
model fitting showed that the best-fit parameters changed as 
follows: AP/foie, respectively Af /back by about ±3 mK, APAfoie 



and AfAback by about ±5°, RMq by about +5 rad m ^. Dfo re- 
mained unchanged. 

It is also interesting to compare the results from our model 
based on the large-scale adjustment of the Effelsberg polariza- 
tion data by the DRAO 1 .41 GHz survey presented in this paper 
with those based on the quite limited Dwingeloo data, which 
were previously used by WoUeben t200U and Wolleben & 
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Reich ( I2004> . For the case of the Taurus region the Dwingeloo- 
derived corrections leaves the EfFelsberg Stokes U data un- 
changed and gives a 138 mK offset at 1.4 GHz for the Stokes 
Q map. When comparing the Dwingeloo and DRAO polar- 
ization intensity scales a difference of 15% was found in the 
sense that the Dwingeloo intensities are too high (Wolleben 
et al., in prep.). In addition the DRAO U and Q maps show 
some gradients or warping across the field due to their better 
sampling compared to the Dwingeloo data. Running the same 
model fitting we get the following differences for the average 
of all clouds analyzed: Af /fore, respectively AP/back, by about 
-25 mK and -12 mk. APAfore and APAback by about 0° or -3°, 
RM() by about -4.3 rad m"^ and DPq - 1 stays unchanged. 

Although it is essential to take the large-scale polarized in- 
tensities for the analysis of Faraday screens into account, their 
uncertainties do not change the qualitative results and the quan- 
titative errors are moderate. The model results for the intrinsic 
RM() of the Faraday screens vary within about 20% and allow 
a conclusive physical interpretation. 

5. Analysis and interpretation 

The parameters for each object we obtained from our model 
fits are shown in Table |2l At a frequency of 1408 MHz, the 
mean values' for the emission component in front of the screen 
is P7foi-e = 175_± 40 mK and PAfoi-e ± 5°. The back- 

ground yields P/back = 142 ± 37 mK and PAback = -1 1 + 9° to 
the observable polarization. The intrinsic rotation measures of 
the Faraday screens average to RMq - (-25 ± 4) radm"^. All 
Faraday screens can be explained with only a small degree of 
depolarization or none at all. 

In the following discussion, the spatial coincidence of po- 
larization minima with the molecular cloud will be interpreted 
as a strong hint that the Faraday effects, which give rise to the 
polarization feature, are associated with the molecular cloud. 
Therefore the distance to the Faraday screens is 140 pc, the 
same as it is for the molecular cloud. 

Our results show that the background component in most 
cases is smaller than the foreground emission. This is an unex- 
pected result since the foreground emission originates within 
only 140 pc. At a Galactic latitude of about -10° the z distance 
is about 24 pc and the emission contributing to the background 
should originate within several kpc. This may indicate that the 
background layer consists of several superposing emission lay- 
ers. The longer the line-of-sight the more likely are changes 
in the Galactic magnetic field direction and/or Faraday rotation 
within the ISM. This means that depolarization effects increase 
with the line-of-sight and the observed polarized intensity de- 
creases. The fit-results seem to reflects this. 

At a distance of about 140 pc of the molecular clouds 
the size of the Faraday screens is of the order of 2 pc, when 
spherical symmetry is assumed. Intrinsic RM() values of up to 
-25 rad mT^ require an excessive thermal electron density or 

' Means were derived by averaging Stokes U and Q separately and 
calculating PI and PA afterwards. The errors of PI and PA were cal- 
culated by the error law of Gauss using the standard deviations of U 
and Q from their means. 



an excessive regular magnetic field component along the line- 
of-sight, when compared to average Galactic values. Since RM 
is the product of both quantities, further constrains are required 
to derive the appropriate numbers. We have checked available 
Ha data from the "full-sky H-alpha map" (Finkbeiner l2003> . 
which do not show any related emission (Icr detection level 
of 0.52 Rayleigh) and constrain to less than 0.8 cm"^ for 
thermal electrons from ionized hydrogen. Upper limits from 
thermal radio emission give < 2 cm With an upper limit 
of 0.8 cm"^ for the thermal electron density n^, a regular mag- 
netic field strength exceeding 20 ;uG within the Faraday screens 
along the line-of-sight is needed. 

Is the existence of such enhanced localized magnetic field 
structures in agreement with the observations? While magnetic 
fields directed towards us will not enhance the total intensity 
synchrotron emission, similar magnetic field enhancements 
should also exist in orthogonal directions. We can roughly cal- 
culate the expected fluctuations of the Galactic background 
emission. Let us assume a magnetic field with an orthogonal 
component of 3 yuG extending for 1 kpc along the line-of-sight 
and giving about 1 K brightness temperature, which is a typi- 
cal value at medium Galactic latitudes at 1.4 GHz. A magnetic 
field enhancement of 20 fiG over 2 pc, keeping the relativistic 
electron density constant, gives about 50 mK for a temperature 
spectral index of yS - -2.7. Such magnetic bubbles will cause 
fluctuations of a few times the confusion limit, which is about 
15 mK Tb for the Effelsberg telescope at 1.4 GHz. Such fluc- 
tuations are quite difficult to identify in practice, and we con- 
clude that the existence of local magnetic field enhancements 
over size scales of a few parsec is in agreement with the ob- 
servations. Bubbles with field strength exceeding about 30 fiG 
or 5 pc in size, however, cause total intensity enhancements of 
the order of 150 mK Tb. Such structures should be clearly de- 
tectable on available radio maps and therefore are unlikely to 
exist in large numbers. 

5.1. The origin of Faraday screens 

Although the origin of the Faraday screens discussed here is 
not clear, we comment on a number of known processes which 
might be able to cause an enhancement of ne and/or the B-field. 

5.1.1. Ionization in a PDR 

In a photodissociation region (PDR) the physical and chemical 
properties of the interstellar medium are dominated by UV ra- 
diation. This is certainly the case for an 0,B star ionizing its 
surroundings, but even the mean interstellar radiation field can 
give rise to PDRs. Most of the atomic and molecular gas in our 
Galaxy is exposed to a far-ultraviolet flux and therefore PDRs 
are expected to be quite common (see HoUenbach & Tielens 
1995). 

Models and observations of PDRs show that with increas- 
ing extinction Ay the chemistry changes from ionized hydro- 
gen at the surface to neutral and finally to molecular hydrogen 
inside the cloud. In a similar sequence carbon exists as in 
the outer region and in form of C and CO inside. According to 
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these models there is a thin layer in the cloud where the ion- 
ization of carbon becomes the main source for electrons, while 
the hydrogen-ionizing photons are already absorbed. 

From Fig. (Sji we roughly estimate an excess reddening 
at the positions of the Faraday screens of Eb-v < 0.5 mag. 
Using the usual gas-to-dust relation A^(Hi + H2)/£'b-v = 
5.8 X 10^' cm^^mag ' (BoWin et al. I1978> . the total neutral 
hydrogen density amounts to n < 470 cm~^ for the size of the 
Faraday screens of 2 pc. 

Several authors have derived a C/H ratio of about 10 
to 4 X 10""*. Variations of C/H associated with the physi- 
cal conditions of the ISM are expected and investigated by 
Gnacinski (2000*). His analysis indicates variations of C/H for 
different fractional abundances of molecular hydrogen /(Hj) = 
2Af(H2)/(A^(H i) + 2NU2- Average values of C/H = 3.55 x lO^^ 
for lines-of-sight with /(Hj) < 10 ^ and C/H = 1.36x10 for 
/(Hj) > lO^^* were found. For extinctions of Ay ~ 3.1 Eb-v < 
1.6 mag, where we observe the Faraday screens, the PDR 
model by Hollenbach et al. a9911 gives /(H,) < 10 I In such 
a case the complete ionization of carbon will enhance the local 
electron density by about 0.2 cm"-'. 

This simple estimate shows that a significant part of the 
thermal electrons in the Faraday screens modeled here may be 
due to carbon ionization in a PDR. In this PDR scenario two 
extreme cases are possible: either no ionization of hydrogen 
takes place and all electrons are from ionized carbon (total ~ 
0.2 cm"^) or both elements contribute to a total electron density 
of «e < 0.8 cm-3(Hii)H-0.2 cmT^C^). With a filling factor of 1, 
a magnetic field strength between 15 //G and 77 fiG is required 
to produce a RMq of -25 radm"-. Further observations of Ho- 
and C^ emission could constrain these limits. 

5.1 .2. Ionization by cosmic rays 

In clouds with moderate gas densities both photoionization and 
cosmic -ray ionization could be important. Theoretical values 
for the resulting ionization are not well known and spatial vari- 
ations are expected. However, referring to McKee (1989 1 a 
molecular cloud is primarily ionized by the external UV radia- 
tion field. Only about a tenth of the mass is ionized by cosmic 
rays. Model calculations by Myers & Khersonsky ( 1995 1 show 
that for gas densities comparable to the ones adopted here, pho- 
toionization at the surface of a molecular cloud exceeds the 
ionization fraction from cosmic rays by a factor of 100. Here 
an exposure to an UV photon flux from the mean interstellar 
radiation field is assumed. Therefore ionization by cosmic rays 
seems to be unimportant. 

5.1.3. Magnetic field enhancement 

A RM enhancement results when the By field gets enhanced. 
There are a number of possibilities to locally enhance the 
regular component of the magnetic field in interstellar space. 
Shock waves from supernova remnants are known to compress 
the field component perpendicular to their velocity direction. 
Stellar winds or expanding Hii regions will also have an ef- 
fect on the ambient interstellar magnetic field. However, there 



is no indication for the presence of a nearby supernova rem- 
nant or a massive hot star with a strong stellar wind in the 
vicinity of the molecular cloud discussed here. Nevertheless 
weak shock waves in the interstellar medium without any clear 
signature may cause some magnetic field compression when 
they hit a massive molecular cloud. Magnetic fields also get 
enhanced during the formation process of molecular clouds, al- 
though missing excessive synchrotron emission limits the mag- 
netic field strength outside the clouds to values not too different 
from those in the interstellar space. Zeeman-splitting observa- 
tions could provide the line-of-sight magnetic field strengths or 
at least an upper limit for the outer regions of molecular clouds. 

6. Local synchrotron emissivity 

The local synchrotron emissivity has been previously deter- 
mined by several authors. Either the low-frequency emission 
in front of Hii regions was measured or model fits of the 
Galactic synchrotron emission were made. The Beuermann et 
al. J1985> model gives about 1 1 K/kpc for 408 MHz, which cor- 
responds to about 0.4 K/kpc at 1408 MHz depending on the lo- 
cal synchrotron spectral index where we assumed /3 = -2.7. The 
Phillipps et al. ( 198 1 ) 408 MHz model of the Galaxy places the 
Sun near a local minimum. They report even lower values cor- 
responding to less than 0.14 K/kpc at 1408 MHz. For the case 
of low-frequency measurements towards optically thick H ii re- 
gions as discussed by Fleishman & Tokarev (1995 1 at 10 MHz 
and Roger et al. ( .1999) at 22 MHz an enhanced local volume 
synchrotron emissivity is indicated. However, this is largely 
based on the nearest H ii region ^-Oph (Sh 27) at ~ 170 pc dis- 
tance, which shows excessive foreground synchrotron emission 
as discussed by Roger et al. ( 1999 1. Excluding ^-Oph, Roger et 
al. d 19991 1 report a mean value of 30 K/pc at 22 MHz, which 
corresponds to 0.9 to 0.6 K/kpc at 1408 MHz for ^ of -2.5 or 
-2.6, respectively. 

Our mean polarized emission at 1408 MHz in front of the 
Taurus molecular cloud is about 0.18 K, which requires a syn- 
chrotron emission of at least 0.24 K for the maximum possi- 
ble percentage polarization. The distance to the Taurus cloud is 
140 pc. Assuming the same distance for the Faraday screens, 
this value corresponds to an emissivity of about 1.7 K/kpc. 
Since this is a lower limit set by the maximum possible 
polarization (~ 75%), the local emissivity is likely signifi- 
cantly higher Other than for the Hii region ^-Oph (Sh 27), 
which is at a similar distance, there is no indication for ex- 
cessive foreground emission towards the Taurus molecular 
clouds. However, the local emissivity may depend on direc- 
tion. This may be investigated by measurements of this kind 
towards molecular clouds in other areas of sky. In case the lo- 
cal synchrotron emissivity derived towards the Taurus cloud 
holds for all directions, synchrotron models of the Milky Way 
(Beuermann et al. 1985 Phillipps et al. 1981) may need mod- 
ifications. Enhanced local emissivity has a significant effect 
when modeling the thick disk properties of the Galaxy. The 
large height of ~5 kpc or more inferred by the Beuermann et 
al. ( I1985> and Phillipps et al. ( I1981> models will be reduced 
or its emissivity decreases. It should be noted that the modeled 
Galactic halo extends several kpc above the plane. Such ha- 
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los are rare and place the Milky Way within a small group of 
galaxies showing similarly huge halos (e.g. Dumke & Krause 

7. Conclusions 

We have used the Effelsberg telescope for observations of the 
polarized emission in the Taurus region and added the miss- 
ing large-scale components by making use of a new DRAO 
1 .4 GHz polarization survey. We observe minima in polarized 
intensity along the boundary of a molecular cloud complex 
located at a distance of 140 pc. Excessive RMs were mea- 
sured based on 1408, 1660 and 1713 MHz data. The polar- 
ization angles do not follow a dependence, which indicates 
the presence of components with different physical proper- 
ties seen in superposition along the line-of-sight. The observed 
systematic variation of polarization angles with polarized in- 
tensity across the Faraday screens were fitted by a model as- 
suming a constant large-scale foreground and also background 
emission. The large-scale components at /i 18 cm were scaled 
from absolutely calibrated All cm data with a spectral index 
of yS^s - -2.7 and /?Mls - radm"^. We calculate physical 
parameters for the magneto-ionic medium at the surface of the 
molecular clouds, which differ largely from average interstellar 
values. Since there is no observational evidence for enhanced 
thermal electron density from available data, an enhanced B- 
field is needed to account for the large RMq observed. Towards 
Faraday screens our model predicts significant variations of 
RMoba for different frequency intervals (see Fig. |5|i. Recently 
published RM maps such as those of Haverkorn et al. (2003 1, 
Gaensler et al. C2001j and Uyaniker et al. (2003) are based on 
fits assuming a A^ dependence of multi-channel data. The fit is 
often imperfect or impossible. This may be due to limitations 
in the measurement accuracy. However, the superposition of 
various components of the magneto-ionic interstellar medium 
results in deviations from a A^ dependence and this should be 
taken into account when quoting RM values. 

From the model fits presented we find quite high values 
for the foreground synchrotron emissivity in the direction of 
the Taurus molecular clouds. If this emissivity represents the 
omnidirectional local emissivity, it is several times larger than 
previously assumed. This will have consequences for the extent 
of the large-scale-halo or the thick disk emission of the Galaxy, 
which will shrink in size or is weaker than previously assumed. 
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